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Abstract. - We demonstrate via several examples that a uniform drift velocity gives rise to 
anomalous aging, characterized by a specific form for the two-time correlation functions, in a 
variety of statistical-mechanical systems far from equilibrium. Our first example concerns the 
oscillatory phase observed recently in a model of competitive learning. Further examples, where 
the proposed theory is exact, include the voter model and the Ohta-Jasnow-Kawasaki theory 
for domain growth in any dimension, and a theory for the smoothing of sandpile surfaces. 



Aging is a characteristic feature of systems such as glasses and spin glasses in their low- 
temperature phase. More generally, aging is a consequence of slow dynamics, which can be 
observed in a broad range of systems far from equilibrium M . Aging phenomena are commonly 
described in terms of the two-time correlation and response functions of the relevant variables, 
which often exhibit scaling behavior in the two-time plane ||. 

In this Letter we demonstrate that this aging scenario is strongly affected by the presence 
of a uniform drift velocity, cither intrinsic to the microscopic dynamics, or imposed by external 
conditions. On general grounds, a uniform velocity field is expected to be relevant for long- 
time dynamics. Consider for simplicity a diffusive relaxational dynamics, characterized by 
a diffusion constant D. To a non-zero velocity V naturally corresponds a finite relaxation 
time, D/V 2 , beyond which the dynamical evolution can be qualitatively different from that 
with V = 0. 

Non-equilibrium phenomena in the presence of velocity fields, either uniform (drift) or 
non- uniform (shear), and more generally in a range of situations with a preferred direction, 
have been described in recent works. For instance, it has been shown that a shear flow induces 
a strong anisotropy in coarsening dynamics |}| . The influence of drift on phase-ordering dyna- 
mics has also been investigated in Ref. The peculiar scaling properties of the flocking 
model of Ref. [|5| are yet another consequence of directionality in non-equilibrium systems. 

Here we show that a variety of models exhibit a novel type of non-equilibrium dynamics 
in the presence of a uniform drift velocity V, characterized by the specific scaling form (|J) for 
the two-time correlation function. We refer to this phenomenon as anomalous aging. Uniform 
drift velocities cannot here be transformed away naively, by invoking Galilean invariance. 
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We first investigate in detail the oscillatory phase described recently in the context of a 
model of competitive learning || . The anomalous aging behavior of this specific model will 
be explained at a quantitative level by the combined effects of drift and diffusion. Anomalous 
aging will then be shown to be a common feature of a range of non-equilibrium models in the 
presence of a drift, including the voter model, the OJK theory for domain growth, and the 
smoothing of sandpile surfaces 

We begin with a brief description of our model of competitive learning Q. Two types 
of individuals, considered as quick and slow learners, are respectively classified by their ef- 
ficiencies (coded as Ising variables r/ x = ±1), and occupy the sites of a regular lattice. The 
updating dynamics of the efficiencies involves two steps. First, a majority rule favors the con- 
version of every individual to the dominant type among its neighbors. In the next step, there 
is a majority rule involving the outcomes of individuals in some stochastic game, favoring 
the conversion of every individual to the type which has had the most favorable outcomes 
among its neighbors at the last iteration. Our model has two variants, the 'interfacial' and 
the 'cooperative', which differ in the precise formulation of the second step. 

We consider the situation where both types of individuals, while distinguishable, are 
equally slow learners, so that no type wins a priori. We refer to this case as coexistence, 
in analogy with ferromagnetic models. Several kinds of ordered and disordered phases, i.e., 
types of non-equilibrium steady states, have been found in the phase diagram of both variants 
of the model [|| ; most phase transitions were found to belong either to the Ising or to the voter 
universality classes. However, an unusual type of coarsening was found in the cooperative ver- 
sion of the model, where the order parameter underwent more or less regular oscillations in 
time, hence the name oscillatory phase. 

As the model of competitive learning is entirely defined by dynamical rules, without refer- 
ence to a Hamiltonian (and therefore without the concept of detailed balance), details of the 
dynamics can influence its steady states. It turns out that the updating scheme, including 
the order in which sites are updated (sequential or parallel, ordered or random), does affect 
the phase diagram. From now on, we focus on the cooperative model at coexistence on the 
square lattice, characterized by a single parameter, < p < 1, which represents an individ- 
ual's probability of obtaining a successful outcome. In Ref. Q ordered sequential dynamics 
were used, and an oscillatory phase was found for p < p Q (p Q « 0.136). On the other hand, 
parallel dynamics leads to a ferromagnetically ordered phase for p < p c (p c « 0.259), while 
random sequential dynamics leads to a disordered phase with very weak correlations between 
efficiencies. 

The special feature of ordered sequential dynamics that is responsible for the oscillatory 
phase turns out to be its directionality. Consider a square sample with sites x = (m, n) with 
m, n = 1, . . . , L. For ordered sequential dynamics, neglecting boundary effects, the next site 
updated after xt = (m, n) is xt+st = (m, n + 1), with St = 1/L 2 . We define the directionality 
of this updating scheme as d = {x. t +st — x t ). We have d = (0,1) for ordered sequential 
dynamics, while d = for parallel dynamics or random sequential dynamics. Other updating 
schemes have been considered, including alternating sequential dynamics: the sample is first 
swept in a regular way, such as e.g. ~x.t+8t = (w, n + 1) if x t = (m, n), and then in a different 
regular way, such as ~x.t+8t = (m,n — 1), or x.t+st = (m + l,n). We have d = in the first 
example, and d = i(l, 1) ^ in the second. The model exhibits an oscillatory phase at low p 
if, and only if, the updating scheme has a net directionality d ^ 0. 

In the following, we restrict the analysis to the prototypical example of ordered sequential 
dynamics. Spin (efficiency) patterns are advected with a constant drift velocity V, due to 
the directionality d of the dynamics. The direction of V is given by the arrow of d, while its 
strength V depends e.g. on the parameter p. If periodic boundary conditions are used, spin 
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Fig. 1 - Logarithmic plot of efficiency correlation function C(t, 0) against t. Full line: least-square fit 
of data for t > 3 (slope = —0.34). 



clusters wrap around the sample, so that the magnetization oscillates more or less regularly, 
with period T\ = 2L/V, which defines a time scale associated with the drift velocity. This 
phenomenon was depicted in Figure 11 of Ref. ||. A correlation analysis of a long magneti- 
zation record allows for the accurate determination of V. We find that V decreases roughly 
linearly with p, from V « 1.14 at p — to V ~ 0.97 at p = p a w 0.136. The drift velocity 
thus remains of order unity throughout the oscillatory phase, including the critical point p , 
where it does not seem to exhibit any kind of critical behavior. 

Next, a coarsening phenomenon takes place if the system starts from a random initial state 
with uncorrelated efficiencies. The commonly accepted picture of phase-ordering dynamics 
with no conservation law [0] involves diffusive behavior, hence the emergence of a second 
time scale, T 2 = L 2 /D, with the diffusion coefficient D being a microscopic constant. The 
competition between these two time scales will turn out to be central to anomalous aging [see 
paragraph below Eq. (||)]. 

Along the lines of investigations of aging properties in glassy systems, we next measured 
the two-time correlation function of the efficiency of a given individual, C(t, s) = (T]x.(s)r] x (t)) , 
with < s (waiting time) < t = s + t (observation time) . We have chosen to work at p = for 
definitencss. For each value of the waiting time s, data were averaged over more than 1,000 
samples of size 300 x 300. The particular case s = corresponds to the correlation of the 
efficiencies at time t with their random initial values. Figures [l] and || respectively show plots 
of C(t, 0) and C(t, s) for s > 0. Figure |l] demonstrates that C(t, 0) decays exponentially, just 
as in any equilibrium process, with a microscopic relaxation time tq w 1/0.34 w 2.9, while 
Figure |^ clearly shows non-stationarity or aging effects: correlations decay more and more 
slowly when the waiting time is increased. 

In order to elucidate these observations, and especially to understand better the role of the 
drift velocity V on the non-equilibrium dynamics of the model, we propose the following phc- 
nomcnological linear theory meant to incorporate the combined effects of drift and diffusion, 
in a natural fashion and at a minimal level of sophistication. The efficiencies ?? x are replaced 
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Fig. 2 - Logarithmic plot of efficiency correlation function C(t, s) against r = t — s, for several values 
of the waiting time s. 



by a continuous field (f>(x,t). The equal-time correlation function C(x, t) — (0(x, i)0(O, t)) is 
assumed to obey the diffusion equation 

dC(x,t)/dt = D 1 W' 2 C(x,t), (1) 

as the drift velocity is expected to have no effect on one-time observables. A disordered initial 
state corresponds to C(x, 0) = (x) at t = 0, in dimension d. The two-time correlation 
function C(x, i, s) = (0(x, i)0(O, s)) is assumed to be affected by both the drift and the 
diffusion mechanisms: 

dC(x, t, s) /dt = D 2 V 2 C{x,t,s) - V • VC(x, *, s) (2) 

(t > s), with initial condition C(x, s, s) = C(x, s) at t = s. The diffusion coefficients D\ and 
Z?2 are possibly different. Convoluting the Green's functions of both equations, and neglecting 
pre-exponential factors, we obtain 

C{t = s + r,s)^^{- A{D ^ D2T) ). (3) 

This analytical result is in excellent agreement with our results above. First, Eq. (||) 
predicts that C(t, 0) decays exponentially, with a relaxation time to = 4D2/V 2 , as seen in 
Figure [l], where the values of V and to yield D2 « 0.94. Second, in order to check the specific 
form of Eq. (||), the data for lnC(t,s) have been plotted as a function of t 2 /(t + as). A 
convincing data collapse and linear behavior are observed (Figure ||) for the optimal value 
a ~ 0.55, yielding D\ = aD 2 ~ 0.52. Moreover, Eq. (||) predicts that the slopes of Figures [l] 
and U are equal; the measured slopes of the least-square fits, —0.33 and —0.34, coincide within 
errors, in agreement with the above. 

Our result (||) is the centerpiece of anomalous aging. The correlation function therein is 
clearly aging, in the general sense that it is not stationary, and decays more and more slowly 
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Fig. 3 - Logarithmic plot of the data of Figure against t 2 /(t + 0.55 s). Full line: least-square fit of 
data with abscissa > 2 (slope = —0.33). 



as the waiting time s is increased. The expression (g) does not, however, exhibit any simple 
scaling form as a function of the time ratio t/s, since it involves the intrinsic time scale 
ro = 4D2/V 2 , emblematic of the competition between drift and diffusion, as expected from 
the discussion in the beginning of this Letter. 

For r <C s, the correlation function in Eq. (||) has a Gaussian scaling form 

C(t,s) « exp ( - t 2 /(t iS )) (r = t - s < s), (4) 

with T\ — AD\/V 2 = tqD\/D2- The correlation is characterized by an effective relaxation 
time 

r cB (s) « (ri S ) 1/2 , (5) 

which diverges, albeit anomalously slowly, with the waiting time s. 

Eq. (H) can be recast as (Vr c g ) 2 w ADis, to be put in perspective with the similar quadratic 
relationship (VT1) 2 = ADT2 = ^L 2 which holds between both time scales introduced above. 
The stationary or time-translationally-invariant (TTI) regime usually observed for r -C s, 
where the correlation function only depends on r, is absent in the present case, having already 
died away before our continuum description applied. The TTI regime is also too small to be 
visible on Figure ||. 

For r 3> s, the result (||) only depends on t, and falls off exponentially, as 

C(t, s) w exp(-r/r ) (r = t-s>s), (6) 

just as C(t, 0). This paradoxical emergence of a translationally-invariant regime for long time 
separations is one of the reasons we refer to this aging as anomalous. The crossover between 
both kinds of behavior (Q) and (^|) takes place for r ~ s, where the correlation is already 
exponentially small, of order exp(— s/tq). 

In agreement with general expectations on aging phenomena the result (0) can be 
recast as C(t,s) w with = exp (2(i/ri) 1 / 2 ) and F(a;) = exp (-(lnx) 2 ). 
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Anomalous aging thus appears as a special case of sub-aging, using the terminology of Ref. , 
as the effective age function h(t) grows much more rapidly than real time t. The effective 
relaxation time (||) is recovered as T e g(t) ~ h(t) / (dh(t) / di) . A similar behavior, with the 
same age function h(t), has been found in the spherical model with random asymmetric 
couplings at zero temperature. The ferromagnetic spherical model with a conserved order 
parameter Jjl"| provides another interesting instance of anomalous aging, involving a scaling 
function similar to the above F{x). 

Since Eq. (||) is a simple consequence of the combined effects of drift and diffusion, it might 
be expected to apply to a variety of situations. Several examples are discussed below, where 
anomalous aging applies exactly in the regime of long times. 

First, our linear theory is exact for the voter model. The dynamics of this model |f2| [l4| 
amounts to saying that a given site, say x, takes the opinion of any of its neighbors, say 
y = x + e, with the lattice vector e being chosen at random with some weight w(e). If all 
vectors e have equal weights (w(e) = 1/z for all e, with z the coordination number), the 
usual voter model is recovered, with its well-known diffusive dynamics [jD^jLjJ. In the general 
(asymmetric) case where some neighbors are preferred over others, patterns of opinions are 
advected with a velocity V = — V e eu>(e). We have shown that the equal-time and two-time 
correlation functions obey Eqs. (yj) and (g) in the continuum limit, with the identity 

D x = 2D 2 , (7) 

so that anomalous aging applies exactly to the voter model. 

Second, the Ohta-Jasnow-Kawasaki (OJK) theory JlSj , a linear approximation commonly 
used in the study of phase-ordering dynamics also yields, in the presence of drift, the 
anomalous aging result The main difference between the OJK theory and our linear 

theory is that linear equations are assumed to be obeyed by the field <p(x, t) itself in the 
OJK theory, and by its correlations in our approach. We have checked that, despite this 
difference, anomalous aging (Eq. (0)) applies exactly to the OJK theory with drift, again with 
the simultaneous validity of Eq. um . 

Our last example relates to the smoothing of sandpile surfaces Q represented in earlier 
work by noisy coupled equations for the fields h(x,t) and p(x,t). Here h(x, t) represents the 
local height fluctuations of the sandpile surface ( caused by clusters of 'stuck' grains) defined 
with respect to a mean surface, while p(x, t) represents the local density of mobile grains 
rolling down the clusters J?]] . The simplest example of these coupled equations is fully linear: 

dh/dt = D h V 2 h + c-Vh + r)(x,t), 

dp/dt = D P V 2 P - c • V/i, 1 ' 

where r/(x, t) is Gaussian white noise, such that (r)(x,t)r](x',t')) = A 2 5<- d \x - x')6{t - f). 
The equation for h(x,t) is nothing but the Edwards- Wilkinson (EW) equation p6] with a 
flow (drift) term added; in the context of sandpiles, this represents the effect of tilt [71. 

In Ref. ||] it was found that while 5Vi(k, t), a one-time quantity, shows only the smoothing 
exponents of the EW equation (i.e., the flow term can indeed be Galilean transformed away in 
this case), S/,(k, u>), a quantity which involves an integration over many times, shows a clear 
crossover from EW exponents to a smoothing state with a = f3 = 0, when the drift velocity is 
large enough. Evaluating the two-time correlation function of the ft,-field, we find 

C h (t,s) = (h(0,t)h(0,s)) = A 2 J^ du (AttDh (t+s- 2u)) ' d/2 cxp (- ^^-fj— \ , (9) 
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In the scaling regime, where both s and t — s = t are large compared to the intrinsic time 
scale To = ADh/c 2 , the integral is dominated by microscopic times u ~ tq. We thus recover 
the anomalous aging result (||), with D\ — 2Dh, D2 = Dh, so that Eq. (fjj) again holds. 

A corollary of this result is that anomalous aging (or, equivalently, the anomalous smooth- 
ing referred to in Ref. ||] ) will manifest itself at experimentally observable times only for large c 
- this, in the sandpile context, would correspond to unstable packings of grains susceptible to 
rolling under tilt . It would be interesting to see if measurements on such a surface indicate 
that it is composed of large domains characterized by EW roughening exponents stacked end 
on end. Physically, this is one possible realization of the anomalous smoothing scenario - 
inter-domain roughness measurements would here show no significant variation in roughness, 
while intra- domain measurements would show EW roughening exponents. 

To sum up, we have demonstrated that a uniform drift velocity is relevant for the long-time 
dynamics of systems far from equilibrium. We have shown that the prediction (^|) of a simple 
linear theory, referred to as anomalous aging, provides an accurate or even exact description 
of the non-equilibrium dynamics of a wide variety of systems, including the oscillatory phase 
observed recently in a model of competitive learning, the voter model, the Ohta-Jasnow- 
Kawasaki theory for domain growth, and a theory for the smoothing of sandpile surfaces. 

# * * 

It is a pleasure for us to thank L. Berthier, L. Cugliandolo, S. Franz, and M. Laessig for 
fruitful discussions. 



REFERENCES 
[1] C.A. Angell, Science 267, 1924 (1995). 

[2] For reviews, see: E. Vincent, J. Hammann, M. Ocio, J. P. Bouchaud, and LF. Cugliandolo, 
in Comp lex Behavior of Gl assy Systems, Springer Lecture Notes in Physics 492, 184 (1997), 



preprint cond-mat / 9607224 ; J. P. Bouchaud, L.F. Cugliandolo, J. Kurchan, and M. Mezard, in 



Spin Glasses and Random Fields, Directions in Co ndensed Matter Ph ysics, vol. 12, ed. A. P. 



Young (World Scientific, Singapore, 1998), preprint :ond-mat/9702070. 
[3] A.J. Bray and A. Cavagna, J. Phys. A 33, L305 (2000); A. Cavagna, A.J. Bray, and R.D.M. 

Travasso, Phys. Rev. E 62, 4702 (2000). 
[4] J. Bragard, P.L. Ramazza, F.T. Arecchi, S. Boccaletti, and L. Kramer, Phys. Rev. E 61, R6045 

(2000). 

[5] O.J. O'Loan and M.R. Evans, J. Phys. A 32, L99 (1999). 

[6] Anita Mehta and J.M. Luck, Phys. Rev. E 60, 5218 (1999). 

[7] Anita Mehta, J.M. Luck, and R.J. Needs, Phys. Rev. E 53, 92 (1996). 

[8] P. Biswas, A. Majumdar, Anita Mehta, and J.K. Bhattacharjee, Physica A 248, 379 (1997). 

[9] A.J. Bray, Adv. Phys. 43, 357 (1994). 
[10] L.F. Cugliandolo, J. Kurchan, P. Le Doussal, and L. Peliti, Phys. Rev. Lett. 78, 350 (1997). 
[11] L. Berthier, Eur. Phys. J. B 17, 689 (2000). 

[12] T.M. Liggett, Interacting Particle Systems (Springer, New- York, 1985). 

[13] J.T. Cox and D. Griffeath, Ann. Probab. 11, 876 (1983); Contemp. Math. 41, 55 (1985); M. 
Bramson, J.T. Cox, and D. Griffeath, Probab. Th. Rel. Fields 77, 613 (1988); J.T. Cox, Ann. 
Probab. 16, 1559 (1988). 

[14] M. Scheucher and H. Spohn, J. Stat. Phys. 53, 279 (1988); P.L. Krapivsky, Phys. Rev. A 45, 
1067 (1992). 

[15] T. Ohta, D. Jasnow, and K. Kawasaki, Phys. Rev. Lett. 49, 1223 (1982). 
[16] S.F. Edwards and D.R. Wilkinson, Proc. Roy. Soc. London A 381, 17 (1982). 



